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PACS 62 . 20 . M — Structural failure of materials 

PACS 63 . 50 . Lm - (Vibrational states in disordered materials) Glasses and amorphous solids 
PACS 63 . 70 . +h - Statistical mechanics of lattice vibrations and displacive phase transitions 

Abstract - At zero temperature and applied stress, an amorphous packing of spheres exhibits, 
as a function of packing fraction, a jamming transition where the system is sensitive to boundary 
conditions even in the thermodynamic limit. Upon further compression, the system should become 
insensitive to boundary conditions but only if it is sufficiently large. Here we explore the linear 
response to a large class of boundary perturbations in 2 and 3 dimensions. We consider each finite 
packing with periodic-boundary conditions as the basis of an infinite square or cubic lattice and 
study properties of vibrational modes at arbitrary wave vector. Our results can be understood in 
terms of competition between plane waves and the anomalous vibrational modes associated with 
the jamming transition; boundary perturbations become irrelevant for systems that are larger 
than a previously identified transverse length that diverges at the jamming transition. 
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Introduction. — At the jamming transition of ideal 
spheres, the removal of a single contact causes the pre- 
viously rigid system to become mechanically unstable 
and crumble [1-3]. Thus at the transition, the replace- 
ment of periodic-boundary conditions with free-boundary 
conditions destroys rigidity even in the thermodynamic 
limit [4-6]. Recognizing that packings at the jam- 
ming threshold are susceptible to boundary conditions, 
Torquato and Stillingcr [7] drew a distinction between col- 
lectively jammed packings, which are stable when the con- 
fining box is not allowed to deform, and strictly jammed 
packings, which are stable to arbitrary perturbations of 
the boundary. Indeed, infinitesimal changes in the shape 
qf the boundary that leave the contact network invariant 
can make a difference; Dagois-Bohy et al. [8] have shown 
that jammed packings with periodic-boundary conditions 
can be linearly unstable to shearing the box. 

At densities greater than the jamming transition there 
are more contacts than the minimum required for stabil- 
ity [1,2]. In this regime one would expect packings that are 
sufficiently large to be stable to changes in the boundary. 
How does the characteristic size for a stable system depend 
on proximity to the jamming transition? The finite-size 
scaling of quantities such as the contact number [3] is gov- 
erned not by the linear system size, L, but by the total 
number of particles in the system, N = L d , for dimension 



d > 2; such scaling is expected for systems at or above 
their upper critical dimension. In contrast, the sensitivity 
to free- versus periodic-boundary conditions is governed 
by a length scale, £*, that diverges at jamming transition. 
For L 3> £* , the system is stable even with free bound- 
aries [4-6]. 

In this paper, we show that stability for a large class of 
boundary perturbations is governed by a separate length 
scale, fr, that also diverges at the jamming transition. 
Packings with L €t are linearly stable with respect 
to these boundary perturbations. We understand this 
as a competition between jamming transition physics at 
low pressures/system sizes, and transverse acoustic wave 
physics at high pressures/system sizes. The two lengths, 
£* and £t, are the same as the longitudinal and trans- 
verse lengths associated with the normal modes of jammed 
sphere packings [9]. 

We analyze systems with periodic boundary conditions 
composed of equal numbers of small and large spheres with 
diameter ratio 1:1.4 all with equal mass, to. The particles 
interact via the repulsive finite-range harmonic pair po- 
tential 

V(n j ) = U 1 - r T-) 2 (!) 



if 



< cry and V(rij) = otherwise. 



Here is the 



distance between particles i and j, Oij is the sum of the 
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particles' radii, and e determines the strength of the po- 
tential. Energies are measured in units of e, distances in 
units of the average particle diameter a, and frequencies in 
units of y^e/ma 2 . We varied the total number of particles 
from TV = 32 to N — 512 at 36 pressures between p = 10 _1 
and p = 10 -8 . Particles are initially placed at random in 
an infinite temperature, T — oo, configuration and are 
then quenched to a T = inherent structure using a com- 
bination of linesearch methods, Ncwtons method and the 
FIRE algorithm [10]. The resulting packing is then com- 
pressed or expanded uniformly in small increments until a 
target pressure, p, is attained. After each increment of p, 
the system is again quenched to T = 0. 

Symmetry-breaking perturbations. The 

boundary conditions can be perturbed in a number of 
ways. The dramatic change from periodic to free bound- 
aries has been studied in refs. [4-6]. Dagois-Bohy et al. [8] 
considered infinitesimal "shear-type" deformations to the 
shape of the periodic box, such as uniaxial compression, 
shear, dilation, etc. Here, we relax the periodic boundary 
conditions by considering a third class of perturbations 
that allow particle displacements that violate periodicity. 
To do this, we treat our system with periodic boundary 
conditions as a tiling of identical copies of the system 
over all space. Thus, an TV-particle packing in a box of 
linear size L with periodic boundary conditions can be 
viewed as the iV-particle unit cell of an infinite hypercubic 
lattice. 

We assume the particles begin in mechanical equilib- 
rium at positions specified by r?, where i indexes parti- 
cles in each cell and /i indexes unit cells, so that r° is 
the equilibrium position of particle i in cell \i. The energy 
of the system to lowest order in particle displacements 
about its minimum value, = r 
be written as. 
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where the sums are over all pairs of particles i/i and jv 
that are in contact. The equations of motion resulting 
from eq. (2) can be solved by a plane-wave ansatz, m^, = 
Re{e.; exp [i(k ■ R^ — u>t)]}. Here e 2 ; is a d/V-dimensional 
polarization vector, k is a d-dimensional wavevector and 
R^ is the d-dimcnsional vector corresponding to the posi- 
tion of cell jj,. This gives the eigenvalue equation, 



X n (k)e ni (k) = Dij(k)e n j(k) 



where, 



A,(fc) = E 



d 2 U 
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(3) 



(4) 



is the dynamical matrix of dimension dN x dN, and n 
labels the eigenvalues and eigenvectors. From eq. (2), the 



frequency of the modes in the nth branch are w„(fc) = 
\/A n (fe) with eigenvector e„i(fe). Note that when k = 
we recover the normal modes for a system with periodic 
boundary conditions that have been widely studied [9,11]. 

With fc allowed to vary over the first Brillouin zone, 
the eigenvectors comprise a complete set of states for the 
entire tiled system. It follows that any displacement of 
particles at the boundary of the unit cell can be written 
as a Fourier series, 



A ki7l e ni (k) exp [i(k ■ R^)} 



(5) 



Therefore, the system will be unstable to some collective 
perturbation of its boundary if and only if there is some 
k and n for which A„(fe) < 0. This procedure allows us to 
characterize boundary perturbations by wavevector. 1 If 
wc find a wavevector whose dynamical matrix has a nega- 
tive eigenvalue, it follows that the system must be unstable 
with respect to the boundary perturbation implied by the 
corresponding eigenvector. 

Our aim is to identify the range of system pressures 
and sizes over which the system is likely to be unstable. 
This amounts to finding the conditions where the lowest 
eigenvalue of the sphere packing is likely to be negative. 
We first consider the so-called "unstressed system," which 
replaces the original system of spheres with an identical 
configuration of particles connected by unstrctched springs 
with stiffness given by the original bonds. We will con- 
struct scaling relations for the lowest eigenvalue of the 
unstressed system as well as for the accompanying shift in 
the eigenvalues upon reintroducing the stress. By finding 
the pressures and system sizes where these two quantities 
are comparable, we will thus determine the scaling of the 
susceptibility of packings to perturbations described by 
eq. (5). 

The unstressed system. — The dynamical matrix 
in eq. (4) is a function of the second derivative of the pair 
potentials of eq. (1) with respect to particle positions. To 
construct the dynamical matrix of the unstressed system, 
we rewrite this as as 



d 2 V dr dr dV d 2 r 
+ 



d 2 V 

dr a drp dr 2 dr a drp dr dr a dr 



(6) 



where V is the potential between particles ifi and jv, 
r = Tj V — and a and /3 are spatial indices. The 
second term is proportional to the negative of the force 
between particles. If this term is neglected, there are no 
repulsive forces and the system will be "unstressed" [5,12]. 
The dynamical matrix of the unstressed system, obtained 
from just the first term in eq. (6), has only non- negative 
eigenvalues because it represents a system of unstretched 



1 Note that shear- type deformations can be considered concur- 
rently with the addition of the term A ■ r?^ to eq. (5), where A is a 
global deformation tensor. This term represents the affine displace- 
ment and is neglected for our purposes. 



p-2 



Stability of jammed packings II: the transverse length scale 




Fig. 1: Dispersion relations along the T — M line for the lowest few branches of 2 different 2-dimensional packings of N = 1024 
particles. The F point is at the Brillouin zone center (fe = 0) and the M point is at the zone corner where the magnitude of 
k is greatest, (a) A "shear-unstable" packing (i.e. a sphere packing that is unstable at low k) at p = 1CF 4 (black) and the 
dispersion relation for the corresponding unstressed system (blue), (b) A "fc > 0-unstable" sphere packing (i.e. a packing that 
is stable near k = but is unstable at higher k) at p = 10 -4 (black) and the dispersion relation for the corresponding unstressed 
system (blue), (c) Comparison between the lowest branch in the sphere packing (black) and its unstressed counterpart (blue) 
for the same system as in (b). The dashed magenta line is the difference between the two eigenvalue branches. 



springs. We will use a subscript "u" (as in e.g., A u ) to refer 
to quantities corresponding to the unstressed system. 

The black curves in fig. 1(a) show the six lowest eigen- 
value branches, A(fc), for a 2-dimensional packing of N = 
1024 disks at a pressure p = 10~ 4 . Also shown (dashed 
blue curves) are the lowest eigenvalue branches for the 
corresponding unstressed system, A u (fc), where the sec- 
ond term in eq. (6) has been omitted. Here, the lowest 
branch of the stressed system has negative curvature at 
k = implying that the system has a negative shear mod- 
ulus. Such shear-type instabilities have been well studied 
previously [8]. In the remainder of this paper, we restrict 
our attention to shear-stable packings that are stable near 
k = but potentially unstable at higher wave vectors. We 
refer to this type of instability as a "fc > instability." 

Figure 1(b) compares the 6 lowest eigenvalue branches 
of a sphere packing (black) with a k > instability to 
those of its unstressed counterpart (dashed blue). The 
lowest branch for the sphere packing has positive curva- 
ture at k = 0, but becomes negative at higher fc. This 
implies that the system is unstable to boundary perturba- 
tions corresponding to eq. (5) over a range of wavevectors. 
In contrast, the unstressed system, which by construction 
must be stable, can have only positive (or possibly zero) 
eigenvalues. In fig. 1(c), the dotted magenta line shows the 
difference between the lowest eigenvalue branch of the un- 
stressed system (blue) and of the sphere packing (black). 

The lowest eigenvalue of the unstressed system. 

— We estimate the eigenvalues of the packing by first 
evaluating the eigenvalues in the unstressed system and 
then considering the effect of the stress term (i.e. the 
second term in eq. (6)). In order to obtain the scaling of 
the k > instabilities, we first estimate the eigenvalues at 
the largest wavevector, k~M, *-e- at the M point located 



at the corner of the Brillouin zone. We then extend the 
argument to the rest of the Brillouin zone. 

For the unstressed system, the lowest eigenvalue at the 
corner of the Brillouin zone can be estimated as follows. 
As fig. 1(b) suggests, the mode structure is fairly straight- 
forward. Low-frequency vibrations arc dominated by two 
distinct classes of modes: plane waves and the so-called 
"anomalous modes" that are characteristic of jammed sys- 
tems [2,4]. The lowest plane- wave branch is transverse 
and parabolic at low fc (see fig. 1(b)): cjt.u ~ ct, u & or 
equivalently 

A T ,u « 4, u fc 2 ~ G u fc 2 , (7) 

where G u is the shear modulus of the unstressed system. 

The eigenvalue of the lowest anomalous mode can be 
understood as follows. Wyart et al. [4, 5] showed that 
the density of vibrational states at fc = for unstressed 
systems, D u (u>), can be approximated by a step function, 
so that D u (uj) ps for u> < u>* while D u (u>) sa const for 
lu > uj*. As suggested by fig. 1, the anomalous modes are 
flat in fe, so this is a reasonable approximation not only 
at fc = but over the entire Brillouin zone. Thus, the 
eigenvalue of the lowest anomalous mode is given by w* 
at any fc. 

Note that if wt,u <C cj*, the lowest branch will maintain 
its transverse-acoustic-wave character and hence will re- 
main parabolic in fc all the way to the zone corner. How- 
ever, when wt, u ^ w*, the lowest mode at the corner 
no longer has plane- wave character because the transverse 
acoustic mode will hybridize with anomalous modes and 
will develop the character of those modes. It follows that 
there is a crossover between jamming physics and plane- 
wave physics when wt,u ~ w* or, equivalently, when 

L « £ T = c T ,uA4. (8) 
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Here £t is the transverse length identified by Silbert et 
al. [9]. 

Near the jamming transition, many properties scale as 
power laws with the excess contact number above iso- 
staticity, AZ = Z — Z^ Q , where Z is the average number 
of contacts per particle and Z?j = 2d(l — 1/N) ps 2d. In 
particular, for the harmonic potentials we consider here, 
a;* ~ AZ, G u ~ AZ, and p ~ AZ 2 for dimensions 
d > 2. (These results are easily generalized to poten- 
tials other than the harmonic interactions used here [2].) 
Equation (8), in combination with cq. (7), predicts that 
the crossover will occur at 



PL 4 



const 



(9) 



for d > 2. 
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Fig. 2: The average eigenvalue of the lowest mode at the cor- 
ner of the Brillouin zone of the unstressed system (blue), as 
well as the average difference between the unstressed system 
and the original packing (magenta), in (a) two and (b) three 
dimensions. The blue and magenta data exhibit collapse as 
predicted by eqs. (10) and (11), with the caveat that we were 
unable to reach the low pressure regime in three dimensions 
where we expect the scaling to be different. Data is only shown 
when at least 20 shear stable configurations were obtained. 

A second crossover occurs at very low pressures and is 
due to finite-size effects that change the scaling of AZ to 
AZ ~ 1/N ~ L~ d in d dimensions, independent of p [3]. 
In this regime, w* and G u remain proportional to AZ, 
and thus also scale as L~ d [3]. We therefore expect the 
lowest eigenvalue of the unstressed system at the corner of 



the Brillouin zone (kni 
regimes. 



y/dn/L) to feature three distinct 



low pressure: A u 
intermediate pressure: A u 
high pressure: A u 



L 



1/N 2 ~ 



-2d 
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(10) 

In two dimensions, we expect that X U L 4 will collapse in 
all three regimes as a function of pL . 

This prediction is verified in fig. 2. The blue symbols 
in fig. 2(a), corresponding to A U L 4 of the unstressed sys- 
tem in d = 2, exhibit a plateau at low pressures/system 
sizes. At intermediate pL A , the blue symbols have a 
slope of 1 and at high pL A a slope of 1/2, as predicted 
by eq. (10). In three dimensions (fig. 2(b)), we observe 
the two higher pressure regimes, with a crossover between 
them that scales with pL 4 , as expected. We did not reach 
the low-pressure plateau regime in three dimensions be- 
cause it is difficult to generate shear stable configurations 
at low pressures. Note, however, we do not expect the 
crossover to the low pressure regime to collapse in d = 3 
with pL A . 

The effect of stress on the lowest eigenvalue. — 

With the behavior of the lowest eigenvalue of the un- 
stressed system in hand, we now turn to the effects of 
stress to explore the behavior of actual sphere packings at 
the zone corner, ku = Vdn / L. The second term in eq. (6) 
shifts the shear modulus to smaller values without affect- 
ing the scaling with pressure, G ~ y/p [13]. It therefore 
follows that G u — G ~ ^fp. In the high-pressure regime, 
where the lowest mode is the transverse plane wave, we 
therefore expect that A u — A ~ G u — G ~ yjp. Likewise, 
the second term in cq. (6) lowers w* 2 by an amount propor- 
tional to p [6, 14, 15], so at intermediate and low pressures 
it follows that A u - A w uj* 2 - lo* 2 ~ p. 

The difference in the lowest eigenvalues of the sphere 
packing and the unstressed system will therefore feature 
two distinct regimes: 



low and int. p: 
high p: 



(A u - A) 
(A u - A) 



P 
P 



i/*L- 



(11) 



The lowest eigenvalue and stability of the origi- 
nal packing. — Comparing these results to eq. (10), we 
see that the lowest eigenvalue of the packing A should be 
positive in the low pressure limit where A u — A -C A u . At 
high and medium pressures, A u and A u — A are comparable 
and obey the same scaling. One would therefore expect 
instabilities to arise in this regime. At high pressures, how- 
ever, we know that A should be positive since A ~ Gk 2 AI 
and shear stability implies G > [8]. Therefore, fluctu- 
ations about the average scaling behavior are most likely 
to drive the system unstable at intermediate pressures. 
Since this regime collapses with pressure and system size 
as pL A , we expect the fraction of systems that are unsta- 
ble to obey this scaling. This prediction is corroborated 
in fig. 3, where we see that the fraction of systems that 
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Fig. 3: The fraction of shear-stable systems in two and three dimensions that are also k > 0-stable. In red we plot the fraction 
of systems that are unstable at the M point while in black we plot the fraction of systems that are stable everywhere. We see 
that both collapse with pL 4 , or equivalently L/£t, with the expected exception of the low pressure regime in three dimensions. 



are stable at the M point (red data) depends on pressure 
p and system size L as expected. 

We have thus far illustrated our reasoning for wavevec- 
tors at a zone corner, which correspond to the smallest 
wavelengths and thus the most drastic perturbations. Sys- 
tems should therefore be more likely to go unstable at the 
zone corner than anywhere else. However, our arguments 
apply equally well to any point in the Brillouin zone. 

This argument is confirmed in fig. 3. To investigate the 
stability over the entire zone, we computed the dispersion 
relation over a mesh in k space and looked for negative 
eigenvalues. Any configuration with a negative eigenvalue 
at any value of k was labeled as unstable. The black data 
in fig. 3 shows the fraction systems that are stable over 
the entire Brillouin zone. It exhibits the same qualitative 
features as at just the zone corner (red data), but with 
fewer stable configurations overall. 

Discussion. — We have shown that shear stable pack- 
ings can be unstable to a class of boundary perturbations 
that correspond to particle vibrations at non-zero k. The 
susceptibility of packings to such perturbations is governed 
by the transverse length scale £t, which corresponds to the 
wavelength of the transverse plane wave with frequency 
uj*; systems are stable when L 3> £t- Also, due to finite- 
size effects, there is a low-pressure regime where systems 
also become stable. 



Note that our scaling arguments for the unstressed sys- 
tem apply equally well to the original packing. However, 
scaling alone cannot tell us the sign of the lowest eigen- 
value. Treating the unstressed and stressed terms of the 
dynamical matrix separately, and exploiting the fact that 
the unstressed matrix is non-negative, is thus necessary 
for understanding the stability of the packing and empha- 
sizes that "fc > 0-instabilities" are caused by stress, not 
the geometry of the packing. 

We can also consider the stability of the system to all 
infinitesimal boundary perturbations. This includes the 
macroscopic shear-type deformations described by Dagois- 
Bohy et al. [8] as well as those studied here. Note that any 
jammed packing prepared with periodic boundary condi- 
tions is stable against infinitesimal compression at all sys- 
tem sizes. Dagois-Bohy et al. [8] found that in two di- 
mensions, the fraction of states that are shear stable col- 
lapses with pL A and approaches 1 at large pL 4 , consistent 
with the speculation that the criterion for shear stabil- 
ity is L ^> £t- Here, we have shown that shear-stable 
packings are stable to boundary perturbations described 
by eq. (5) for L i?. Moreover, the value of the lowest 
eigenvalue exhibits scaling collapse with L/£t, except at 
very low pressures where such scaling is expected to break 
down due to finite-size effects. 

The combination of these two results suggests that the 
system should be stable to all infinitesimal boundary per- 
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turbations for L ^> £t- This implies that the closer the 
system is to the jamming transition the larger the sys- 
tem must be in order to be insensitive to changes in the 
boundaries. It also implies that the distinction between 
collectively and strictly jammed [7] packings disappears 
when L 3> £t and thus vanishes in the thermodynamic 
limit for any nonzero pressure. 

While the sensitivity of jammed packings to infinitesi- 
mal changes to the boundary is controlled by the diverging 
length scale £t, Wyart el al. [6] argued that the stability 
to a more drastic change of boundary conditions, in which 
the periodic boundaries are replaced with free ones, is gov- 
erned by the larger length £* . Jammed packings with free 
boundaries are stable only for L 3> l\, [6]. Goodrich et 
al. [6] have shown that this length is equivalent, not only 
in scaling behavior but also in physical meaning, to the 
longitudinal length £^ = uj* jcu that was proposed along- 
side £t by Silbert et al. [9] and is the wavelength of the 
longitudinal plane wave with frequency ui*. 

The two length scales £* = £\, and £t represent the sta- 
bility of jammed packings to different classes of boundary 
changes and are both related to the competition between 
plane- wave physics and jamming physics, in that both are 
related to the characteristic frequency w* via the appropri- 
ate sound speed. The two lengths can therefore be viewed 
as central to the theory of the jamming transition. 

To understand why some boundary perturbations are 
controlled by £t while others are controlled by £lj we must 
understand the nature of the transverse and longitudinal 
deformations allowed. In the case of systems with free 
boundaries, the system size of £l needed to support com- 
pression is also sufficient to support shear since £\, > £t- 
The minimum system size needed to support both is there- 
fore £jj. In the case of periodic boundaries, we have shown 
here that jammed systems are stable to transverse pertur- 
bations for sizes larger than £t>. However, they are stable 
to longitudinal perturbations at any system size. The min- 
imum system size needed to support both is therefore £t. 

The same reasoning can be applied to understanding 
the stability with respect to a change from periodic to 
fixed boundary conditions. In that case, neither longitu- 
dinal nor transverse deformations are allowed and it fol- 
lows that systems of any size will be stable to this change. 
Mailman and Chakraborty [16] have studied systems with 
fixed boundary conditions to calculate point-to-set corre- 
lations. Their analysis, as well as arguments based on the 
entropy of mechanically stable packings, reveal a correla- 
tion length that scales like I*. 

The fact that the diverging length scales control the 
response to boundary changes but do not enter into the 
finite-size scaling of quantities such as the contact num- 
ber and shear modulus [3] is consistent with the behavior 
of a system that is at or above its upper critical dimen- 
sion. The fact that power-law exponents do not depend on 
dimensionality for d > 2 is also consistent with this inter- 
pretation [1, 17]. However, critical systems are generally 
controlled by only one diverging length scale. Jamming 



is thus a rare example of a phase transition that displays 
two equally important diverging length scales. It remains 
to be seen whether other diverging lengths that have been 
reported near the jamming transition [1, 16, 18-20] can also 
be associated with boundary effects. 

We thank David L. Johnson for prompting us to 
study vibrations of jammed tilings, Justin Burton for 
suggestions regarding energy minimization methods, as 
well as Tom Lubensky and Anton Souslov for stimu- 
lating discussions. This research was supported by the 
U.S. Department of Energy, Office of Basic Energy Sci- 
ences, Division of Materials Sciences and Engineering un- 
der Award DE-FG02-05ER46199 (AJL,SS,OK) and DE- 
FG02-03ER46088 (SN,OK), and by the NSF via the grad- 
uate research fellowship program (CPG) and the UPenn 
MRSEC DMR11-20901 (SS). 

REFERENCES 



[1] 
[21 
[3 
[4 
[5 
[6 
l~ 
[8] 

[9] 
[10 

[11 

[12 
[13 

[14 
[15 

[16 

[17 

[18 

[19; 
[20; 



O'Hern C. S., Silbert L. E., Liu A. J. and Nagel 
S. R., Phys. Rev. E., 68 (2003) 011306. 
Liu A. J. and Nagel S. R., Annu. Rev. Condens. Matter 
Phys., 1 (2010) 347. 

Goodrich C. P., Liu A. J. and Nagel S. R., Phys. Rev. 
Lett, 109 (2012) 095704. 

Wyart M., Nagel S. R. and Witten T. A., Europhys. 
Lett, 72 (3) (2005) 486. 

Wyart M., Silbert L. E., Nagel S. R. and Witten 
T. A., Phys. Rev. E., 72 (2005) 051306. 
Goodrich C. P., Ellenbroek W. G. and Liu A. J., 
arXiv:1301.6981 [cond-mat.soft]. 

Torquato S. and Stillinger F. H., J. Phys. Chem. B, 
105 (2001) 11849. 

Dagois-Bohy S., Tighe B. P., Simon J., Henkes S. 
and VAN Hecke M., Phys. Rev. Lett, 109 (2012) 095703. 
Silbert L. E., Liu A. J. and Nagel S. R., Phys. Rev. 
Lett, 95 (2005) 098301. 

Bitzek E., Koskinen P., Gahler F., Moseler M. and 
Gumbsch P., Phys. Rev. Lett, 97 (2006) 170201. 
Silbert L. E., Liu A. J. and Nagel S. R., Phys. Rev. 
E, 79 (2009) 021308. 

Alexander S., Physics Reports, 296 (1998) 65. 
Ellenbroek W. G., Zeravcic Z., van Saarloos W. 
and van Hecke M., Europhys. Lett, 87 (2009) 34004. 
Wyart, M., Ann. Phys. Fr., 30 (2005) 1. 
Xu N., Wyart M., Liu A. J. and Nagel S. R., Phys. 
Rev. Lett, 98 (2007) 175502. 

Mailman M. and Chakraborty B., J. Stat. Mech., 
2011 (2011) L07002. 

Parisi G. and Zamponi F., Rev. Mod. Phys., 82 (2010) 
789. 

Drocco J. A., Hastings M. B., Reichhardt C. J. O. 
and Reichhardt C, Phys. Rev. Lett, 95 (2005) 088001. 
Olsson P. and Teitel S., Phys. Rev. Lett, 99 (2007) 
178001. 

Vagberg D., Valdez-Balderas D., Moore M., Ols- 
son P. and Teitel S., Phys. Rev. E, 83 (2011) 030303. 



p-6 



